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ABSTRACT 

We determine the spectrum of particles accelerated at shocks with arbitrary 
speed and arbitrary scattering properties for different choices of the equation of 
state of the downstream plasma. More specifically we consider the effect of energy 
exchange between the electron and proton thermal components downstream, and 
the effect of generation of a turbulent magnetic field in the downstream plasma. 
The slope of the spectrum turns out to be appreciably affected by all these 
phenomena, especially in the Newtonian and trans-relativistic regime, while in 
the ultra-relativistic limit the universal spectrum s ~ 4.3 seems to be a very solid 
prediction. 

Subject headings: cosmic rays - shock waves - relativistic shock jump conditions 



1. Introduction 



The mechanism of diffusive particle acceleration at shock fronts is one of the best studied 
acceleration processes possibly at work in a wide variety of astrophysical environments. After 
the pioneering work carried out in the 70's and early 8 's on non-r elativistic shock fronts 
fcrvmskii 19771 : iBlandford fc Ostriker 19781 : iBell 1978al : iBell 1978bh and sl ightly later on 



the extension of the ac celeration process to the case of relativistic shocks (jPeacock 1981 



Heavens &: Drury 19881 ) numerous developments took place mainly driven by three scientific 
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goals: 1) relax the assumption of test particles, therefore introducing the possibility of shock 
modification due to the dynamical reaction of the accelerated particles on the shock itself (see 
(IMalkov fc Drury 200ll ) for a recent review); 2) introduction of self-generation of scattering 
as a result of streaming instability; 3) set the bases for a general theory that could be applied 
to shocks with arbitrary speed. 

From the phenomenological point of view these efforts are motivated by the wealth of 
situations in which the acceleration process can potentially be at work, from the case of 
supernova remnants (SNRs) to that of gamma ray bursts and more in general to the case 
of relativistic jetted plasma ejections from compact objects. The case of SNRs is especially 
relevant for its potential implications for the origin of galactic cosmic rays. 

The point 1) above has received much attention: the problem o f dynamical reaction 
of the accelerated p articles has bee n faced by using two-fluid models ( Drury &: Volk 198ol; 
Drurv fc Volk 198lh. kinetic models jMalkov 199 71: IMalkov et al. 200ol : Blasi 20021 ; blasi 2004 
Blasi et al. 20051 : Amato fc Blasi 2005h and numerical approaches, both Monte Carlo and 
other simulation procedures (jJones fc Ellison 199ll ; iBell 19871 : lEllison et al. 19901 : lEUison et al. 1995 



Ellison et al. 19961 : IrCang fc Jones 19971 : IrCang fc Jones 20051 : IrCang et al. 2002h . All these 



approaches are however limited to non-relativistic shocks. The only case that we are aware 
of in which shock modification has been calculated for a relativistic shock is the numerical 



simulation of (lEUison &: Double 2002l ). Recently (lAmato &: Blasi 20061 ) have presented an 
semi-analytical solution of the combined problem of particle acceleration at modified shocks 
with self-generated Alfven waves (point 2), though treating the wave amplification in the 
context of quasi-linear theory. 

The point 3) above, namely the construction of a unified description of the accelera- 
tion process applicable to any shock, has been successfully solved, in the framework of a 
test particl e approach, by (Vietri 20031; Blasi fc Vietri 2005 ). and applied to some cases of 
interest by iMorlino et al. 2006h . falasi fc Vietri 2005h illustrated a wide variety of compar- 
isons of the theory with existing work in both cases of non-relativistic and ultra-relativistic 
shocks. In t he latter case, much work has been done by several authors, who used numeri- 
cal methods (IBednarz fc Ostrowski 19981: lAchterberg et al. 200ll : iLemoine fc Pelletier 20031: 

Niemiec fc Ostrowski 200*41" Lemoine fc Revenu 2006 ) or semi- analytical approaches ( Kirk &: Schneider 198 
Gallant fc Achterberg 19991 : iKirk et al. 2~00~oL though limited to specific situations (for in- 
stance small pitch angle scattering). 



This paper is t he third of a series, in which the original calculations of (jVietri 2003 



Blasi fc Vietri 20051 ) are applied to investigate cases that can be of some phenomenological 
relevance. Here we determine the changes to the spectral index of the accelerated parti- 
cles, taking into account the possibility that different physical situations may change the 
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equation of state of the plasma downstream of the shock. In particular collective plasma 
effects in the proximity of a collisionless shock may be responsible for different levels of ther- 
malization of the electron and proton gas components behind the shock. These effects are 
expected to provide a channel of energy transfer from protons to elec trons (see for example 



(IBegelman & Chiueh 19881 : iHoshino et al. 1992 ; iGallant et al. 19921 )) though not through 



particle-particle collisions. Other collective processes might result in the generation of a 
downstream turbolent magnetic field, which in turn may change the equation of state of the 
downstream fluid. We describe these effects in a phenomenological way, by introducing a 
parametrization of the energy exchange between electrons, protons and magnetic field energy 
densities in the downstream plasma, and we use such parametrization to infer the changes 
in the equation of state and the spectral slope of the particles accelerated at the shock front. 

The paper is organized as follows: in §2] we describe the jump conditions for a strong 
shock moving with arbitrary speed. In §3]we obtain the equation of state for the downstream 
plasma for different cases of interest. We present our results on the spectrum of accelerated 
particles in ^ and conclude in S|5l 



2. Relativistic jump conditions for strong shocks 

The jump conditions describing the conservation of mass, momentum and energy across 
a shock front moving with velocity V s h = 0sh.c in a medium with density ni, pressure p± and 



energy density e\ are as follows {e.g. (IKirk fc Duffy 19991 ; IGallant 20021 )): 



Ti/3ini = T 2 (3 2 n 2 (1) 

r?A( ei + Pl ) = rjp 2 (e 2 + P2 ) (2) 

TlPl{e 1 + Pl )+p 1 = T 2 2 Pl(e 2 +p 2 )+p 2 . (3) 

Number densities (n), pressures (p) and energy densities (e) are all measured in the comoving 
frame of the plasma they refer to, while the Lorentz factor I\ (T 2 ) of the upstream (down- 
stream) plasma are measured in the shock frame (clearly Ti = T s h = (1 — /5^)~ 1//2 ). The 
indexes '1' and '2' refer to the upstream and downstream plasmas respectively. 

The equations above can be easily generalized to the case of presence of non negligible 
magnetic fields upstream, but in the following we shall assume that the dynamical role of such 
fields is always fully negligible, and shall therefore ignore the corrections in the conservation 
equations. 

The system of equations Eqs. ([I])-© can be solved once an equation of state for the 
plasma has been fixed, in the very general form e = e(n,p). For simplicity, in the following 
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we shall limit ourselves with the case of strong shocks, namely shock waves propagating in 
cold pressureless media, so that p\ = and e\ ~ r^mc 2 = pic 2 . In this framework Eqs. (j2J) 
and ([3]) become: 



r^ lPlC 2 = T 2 2 f3 2 (e 2 +p 2 ) and (4) 
r 2 /5 2 Pl c 2 = T 2 2 (3 2 2 (e 2 +p 2 )+p 2 . (5) 

We assume that the equation of state has the form e 2 = p 2 c 2 F(p 2 / p 2 c 2 ), or in terms of 
the normalized variables p 2 = and e 2 = e 2 / p 2 c 2 : 



h = F{p 2 ). (6) 

Most cases of astrophysical interest are well described by this functional form for the 
equation of state of the downstream g discussed in Sec. El 

Using the equation of mass conservation, Eq. (pQ), Eq. (j4j) becomes 



r 2 = ?i/(e 2 {p 2 ) +p 2 ) = gi(p 2 ), (7) 

while from ([5]) we have 

= (4 - i)/(e 2 fe) 2 - Pi - 1) = g 2 2 (p2) ■ (8) 

The solution for p 2 can be obtained solving numerically the equation gi(p 2 ) = fi^fe)- Once 
p 2 is known, the equation of state gives e 2 while Eq. ((7j), or equivalently Eq. ([8]), gives T 2 . 
Finally Eq. (CQ) gives the number density n 2 . At this point it is also easy to determine the 
velocity jump and the shock rg = /3i//3 2 , which is a crucial parameter for the description of 
the process of particle acceleration at the shock front. 



3. Equations of state for the downstream plasma 

In this section we consider several instances of equations of state for the downstr eam 
gas, in addition to the well known and widely used Synge equation of state (jSynge 19571 ). In 
§3.11 we discuss the case of a downstream plasma made of two independent particle species 
that may thermalize to different temperatures. In §3.21 we introduce the possibility that the 
the proton and electron components are coupled in a collisionless way. In §3.31 we discuss 
the modification of the equation of state due to generation of a turbulent magnetic field in 
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the downstream plasma. Finally in §3.41 we make an attempt to consider the most general 
case in which all the effects described above are taken into account. 



3.1. The case of a plasma with independent particle species 



An equation of state which is widely used in the literature was introduced by (jSynge 19571 ) 
The basic assumption is that the plasma consists of a single component with temperature 
T, and 

e + p = pc 2 G (mc 2 /k B T) (9) 

where G(x) = K 3 (x) / K 2 (x) and K 2 , K 3 are the modified Bessel functions. For x S> 1 {i.e. 
kv,T <C mc 2 ) Eq. (jUJ) reduces to the classical Newtonian equation of state, e = pmc 2 + 3p/2, 
while in the opposite limit i <C 1 (i.e. k^T ^> mc 2 ) the ultra-relativistic equation of state 
is recovered, e = 3p. 

If the downstream plasma can be well described as an ideal gas (p = nksT) made of a 
single component, Eq. (jSj) can be rewritten in terms of the normalized variables introduced 
in the previous section: 

e = G(l/p)-p, (10) 

which has the functional form assumed in §2]with F(p) = G (1/p) —p. The Synge equation of 
state describes correctly the behaviour of the plasma in the ultra-relativistic and newtonian 
limits: applying the procedure illustrated in §2] we easily find the velocity jump rp, which 
is plotted as a solid line in Fig. [TJ as a function of the product T s h/3 s h- In the case of 
strong Newtonian shocks the well known result rp = 4 is recovered. In the limit of an ultra- 
relativistic shock, when the downstream fluid obeys the ultra-relativistic equation of state 
for the gas, rp = 3. 

The generalization of the Synge equation of state to the case of two (or more) indepen- 
dent particle species with temperatures Tj is rather straightforward. Of particular interest is 
the case in which the temperature of the i—th species is simply due to the isotropization of 
the velocity vectors at the shock surface. In this case the energy density of the i—th species 
in the downstream plasma can be written as 62a) = T re \nimiC 2 , or, in terms of dimensionless 
variables: 

£2(2) = r re i, (11) 

independent of the type of particles. Since the normalized energy density is the same for all 
species, the normalized pressures need to be the same too. It follows immediately that for 



- 6- 



the system as a whole one can write: 

e T +Pt = ^Y jPl c 2 G{m l c 2 /k B T l ) = G(l/p T ) , (12) 
PtC — 

i 

where the total quantities have the subscript t- 

3.2. Coupling between thermal protons and thermal electrons 

The formation of collisionless shocks, both in the relativistic and newtonian regime still 
represents a subject of active investigation, in that the mechanisms that allow for an efficient 
transport of information among the particles in the plasma through the exchange of MHD 
waves are poorly known. On the other hand we know that such shocks do exist, which can be 
interpreted as an indirect proof of the importance of collective effects in collisionless plasmas. 
The same type of effects may also be responsible for total or partial thermalization of the 
different components of a plasma. Whether electrons and protons downstream of the shock 
front are in thermal equilibrium or not is a matter of debate. Most likely the answer depends 
on the specific conditions behind the shock of interest. In addition to the thermalization of 
the species, there are several other problems related to our ignorance of the complex physics 
that rules these effects: for instance even the spectrum of the thermal distribution of protons 
(and respectively of electrons) may not be a typical Maxwellian, in particular if a background 
magnetic field makes the distribution of energy in the waves anisotropic. In these conditions 
one should probably introduce a plasma temperature along the field and perpendicular to it. 
The problem of the thermalization of the plasmas around collisionless shocks is also related 
to the issue of the thickness of a collisionless shock, which is usually assumed to be of the 
order of the gyration radius of the thermal proton component. The same collective effects 
also determine the efficiency of injection of particles in the acceleration cycle: it appears 
intuitively clear that in a collisionless shock the processes of thermalization and particle 
acceleration to non-thermal energies are intimately related to each other. 

Lacking a true theory of collisionless energy transport, we adopt here a phenomeno logi- 
cal approach in that we parametrize the degree of equilibration between electron and proton 
temperatures in the downstream region by introducing a parameter £ e , such that the tem- 
peratures of electrons and protons satisfy the relation T e = £ e T p . On very general grounds 
we expect £ e < 1, at least close to the shock, before any collisional effects may possibly 
equilibrate the two temperatures where the plasma has moved away from the shock front. 

If we assume that the electron and proton gas separately behave as perfect fluids, the 
pressures of the two components are related through p e = nk^T e = (, e p P - The equations of 
state of electrons and protons are easily found to be 
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(13) 



(14) 



Here all the normalized quantities refer to the total matter density, namely x = x/n(m e + 
m p )c 2 . In order to solve the equations for the jump conditions at the shock surface we also 
need the total normalized energy densities and pressures: 



We found the appropriate solutions for the jump conditions for several values of the 
parameter £ e in the range m e /m p < £ e < 1. Fig. [T] shows the velocity ratio rp as a function 
of the product T sh /3 sh for £ e = m e /m p , 0.1, 0.3, 0.5, 1 and oo. The case considered in 
section [37T1 corresponds to £ e = m e /m p , while £ e = 1 represents the limit case with protons 
and electrons thermalized at the same temperature. 

The most peculiar feature of the curves plotted in Fig. [T] is the presence of a peak 
approximately at the place where the shock becomes trans-relativistic. The nature of this 
peak and the phenomenological implications of its presence are rather interesting. In order to 
understand the origin of the peak we consider the unphysical case £ e = oo, which corresponds 
to completely cold ions in the downstream plasma. For all curves in Fig. [I] we can clearly 
identify three regimes: 1) both electrons and protons downstream are non-relativistic; 2) 
electrons are relativistic while protons are still non-relativistic; 3) both electrons and protons 
are relativistic. 

The case £ e = m e /m p (or equivalently £ e ~ 0) is the case that is usually studied in the 
literature. In the limit of non relativistic strong shocks this case leads to compression factor 
that asymptotically approaches 4. The derivative of the compression factor with respect to 
(3 s h in this non relativistic regime is zero, as can be shown by using the Taub conditions Eqs. 
CO)-© (or equivalently, and more easily, the non relativistic version, the Rankine-Hugoniot 
jump conditions). On the other hand, if one expands the function G(x) in the equation of 
state of electrons to second order in the variable 1/x (i. e. for mildly relativistic temperature), 
keeping the protons non relativistic, the resulting total equation of state reads 




(15) 
(16) 





15 e e 




(17) 
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It is easy to show that the Rankine-Hugoniot relations give now a compression factor with a 
positive derivative with respect to (3 s h, for small values of [3 s h. Moreover the asymptotic value 
of rp for ultra- relativistic shocks is always 3. This implies that at some point in between 
the compression factor has a peak where the derivative is zero and the compression factor is 
maximum. This is clearly seen in Fig. [TJ in the trans-relativistic regime, electrons become 
relativistic downstream before protons do (and even if the shock is not fully relativistic) 
thereby making the downstream fluid more compressive. To show that this is the correct 
interpretation, one can estimates the value of rg at the peak using the equation of state for 
non relativistic protons (e p = p p c 2 + 3p p /2) and fully relativistic electrons (e e = 3p e ). Under 
this assumption the total normalized energy reads 

_ _ 3 l + 2£ e _ 

Inserting this equation in the Rankine-Hugoniot relations one obtains the following expres- 
sion for the compression ratio: 

" = T+!7' (19) 

Substituting £ e with the values listed in Fig. (JTJ, one recovers values of the compression 
factor close to that of the peaks within an error < 4%. 



3.3. Turbulent magnetic field production 

It is interesting to investigate the possibility that part of the ram pressure of the up- 
stream fluid may be converted into a turbulent magnetic field in the downstream region. 
Here we investigate this scenario and in particular we calculate the compression factor at 
the shock and the spectrum of the accelerated particles. In order to take into account the 
dynamical effect of the turbulent magnetic field it is necessary to generalize the Taub con- 
ditions at the shock, (EJ-Q, by introducing the proper components of the electromagnetic 
stress tensor T^ v . The specific energy density in the form of turbulent magnetic field is 
e m = T 00 while the pressure in the direction identified by the index % is p m j = T n . 

In the plasma reference frame, where no electric field is present, the electro-magnetic 
energy tensor can be written as follows: 

= i- {F» a F v a - ri^F a(i F a p) = (b^B u - ^B 2 rTj • (20) 

Here B^ = (0, B) and B 2 = B 2 X + B 2 y + B 2 Z . To be included in Eqs. (TJ])-© has to 
be expressed in the shock frame, where both the scaler B and the field component along 
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the propagation direction B x remain unchanged. Hence the energy density and the pressure 
along the shock propagation direction are 

e m = (B 2 x + B 2 y + B 2 z )/8n, (21) 
p miX = (-B 2 x + B 2 y + B 2 z )/87T. (22) 

Two situations may be of interest here: 1) the turbulent field is created directly behind 
the shock. In this case the strength of the field is equally distributed among the three spatial 
dimensions. 2) the turbulent field downstream results from the compression of a turbulent 
field upstream. In this second case the two components of the field which are perpendicular 
to the shock normal are amplified at the shock while the parallel component is left unaltered. 
In the former case the relation between energy density and pressure is easily obtained to be 

Pm = £ m /3, (23) 

where the factor 1/3 suggests that the magnetic field behaves like a relativistic gas irre- 
spective of the shock speed. On the other hand, in the latter case, if the shock is ultra- 
relativistic, the parallel component of the turbulent magnetic field is negligible with respect 
to the perpendicular components due to the shock compression. If the parallel component is 
neglected, the relation between energy and pressure of the magnetic field is easily obtained 
to be p m = e m . 

In the following we limit ourselves with the case p m = e m /3, but we discuss the case 
p m = e m in Sec. U 

If £ m is the fraction of magnetic energy density with respect to the proton kinetic energy 
density, we can write 

e m = i m (e p - pp) . (24) 

When the magnetic energy equals the kinetic energy of protons, the magnetic pressure 
is smaller than that of protons in the Newtonian limit. On the other hand the two pressures 
become equal in the ultra- relativistic limit. It follows that we expect the compression ratio to 
increase in the Newtonian limit as the magnetic contribution increases, while the compression 
factor levels off when the relativistic regime is approached. The total downstream equation 
of state when only protons and a turbulent magnetic field are taken into account is: 

£t = e P + e m = e P (l + £m) - imPp ■ (25) 

Fig. [2] shows the compression factor rp for this situation, for several values of £ m , in the 
range < £ m < 1. The velocity compression factor ranges from 5.0 in the Newtonian limit 
to 3.0 in the ultra-relativistic limit, when protons and the magnetic field are considered in 
equipartition (i.e. £ m = 1). 



2 t i i I i i i i i i i i 1 

0.01 0.1 1 10 100 

^shock Pshock 



Fig. 1. — Velocity compression factor at the shock when the downstream protons transfer a 
fraction £ e of their thermal energy to electrons. 





1 

^shock Pshock 



Fig. 2. — Velocity compression factor at the shock when a turbulent magnetic field is present 
inside the electron-proton plasma, and £ m is the fraction of magnetic field energy density 
with respect to the protons kinetic energy. 
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3.4. The general case 

As a natural conclusion of our exercise, we consider now the case of an electron-proton 
plasma with all the effects introduced above. The total energy of the system is given by the 
sum of the Eqs. f|T3|) . ffT4"|) and ff24l) . The total pressure is the sum of the three contributions 
due to protons, electrons and magnetic pressure. We can write: 



e p + e e + e T 



e P (l + 6 m) + e e - p P ir, 



[l+im)PpG 



pp 
pp 



+ PeG 



Pt(p p 



-Pp(l + 6m + 6e) - Pp&, 
Pp + Pe + Pm 



Pp 



1 i /- £m \ . 6m - 
l+6e- y + T P P 



iePp 



p v 



- i 



(26) 



(27) 



where we also introduced the proton and electron normalized densities p p = m p j (m p + m e ) 
and p e = m e / (m p + m e ). 

In order to avoid the effects of proliferation of free parameters, in the following we limit 
ourselves with a sort of equipartition situation, in which £ e = 6m = 6- Fig- El shows the 
resulting compression factor rp for m e /m p < 6 < 1. 

It is important to stress that for shocks in the newtonian and trans-relativistic regime, 
the magnetic field and the thermal electrons both result in making the plasma more com- 
pressible (the compression factor is as high as 5.6 at the peak T sh (3 sh ~ 0.3, when equipar- 
tition 6 = 1 is assumed). On the other hand for highly relativistic shocks all the three 
components behave in the same way resulting in a compression factor equal to 3. The elec- 
tron contribution turns out to be especially important in the intermediate velocity range 

(o.i<r s /3 s <i). 



4. The particles' spectrum 

The spectrum and angular distribution of the particles accelerated at shocks with ar- 
bitrary speed and arbitrary scattering properties of the fluid can be calculated following 
the theory of particle acceleration put forward by (IVietri 20031 ). This approach requires the 
calculation of the conditional probabilities of a particle returning from upstream or down- 
stream at some direction given the entrance direction. These conditional probabilities were 
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described by (IBlasi &: Vietri 20051 ) in terms of two non-linear integral equations, that were 



solved iteratively 



The theoretical approach of (jVietri 20031 ) and (IBlasi fc Vietri 20051 ) applies equally well 



to cases of Small Pitch Angle Scattering (SPAS) and Large Angle Scatte ring (LAS), and 



to the case of a large scale coherent field upstream (IMorlino et al. 20061 ). Here we limit 
ourselves with considering only two situations, namely that of SPAS, both upstream and 
downstream, and that of a large scale fi eld upstream, with o rientation perpendicular to the 



shock normal (perpendicular shock). In (IMorlino et al. 20061 ) it was shown that the spectral 



shape does not change dramatically with the inclination, with the exception of the cases 
where the shock is quasi-parallel, but these cases lead to insignificant acceleration and are 
therefore physically irrelevant. 

Before showing our results for the different equations of state discussed above, it is 
useful to show the spectral slope and the distribution function of the accelerated particles 
for different values of the velocity compression factor. In Fig. H] we show the spectral slope 
as a function of the compression factor for several values of T s hf5 s h, ranging from 0.05 to 5.0. 
The solid line corresponds to the case of a strong newtonian shock, s(rp) = (rg + 2)/(rg — 1). 
In Fig. [5] we also show the distribution function of the accelerated particles as a function of 
the direction fi measured in the downstream frame. 

As discussed in § §3.1H3.4[ the basic effect of changing the equation of state is to change 
the compression factor at the shock and thereby the shape of the spectrum of the accelerated 
particles. In the following we discuss separately the effect of the interaction between electrons 
and protons in the downstream plasma, and the effect of the generation of magnetic field. 
Finally we shall use Eqs. (|26|) and (}2~Tj) in order to quantify the combined effect of the two 
phenomena. 

Fig. [6] shows the spectral slope as a function of the shock speed for different values 
of the parameter £ e which characterizes the degree of coupling between thermal electrons 
and thermal protons, as introduced in § §3.11 The left panel refers to the case of SPAS both 
upstream and downstream, while the right panel refers to the case of a regular perpendicular 
field upstream (and SPAS downstream). 

Some comments are in order: for very non-relativistic shocks (not shown in the plot) 
the spectrum has the standard slope s = 4. In the ultra-relativistic limit universality is also 
reached, being the results independent on the value of £ e . Most differences in the slope of the 
spectrum of accelerated particles is present for trans-relativistic shocks: a minimum in the 
slope appears for these shocks, deeper for larger values of £ e . This flattening of the spectrum 
is due to the increasing compression ratio at trans-relativistic shock speed, as explained in 
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0.01 0.1 1 10 100 

^shock Pshock 



Fig. 3. — Velocity compression factor when the equation of state of the downstream plasma 
takes into account the exchange of energy between electrons and protons, and the generation 
of a turbulent magnetic field. Here we assumed £ e = £ m = £. 
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Fig. 4. — Slope of the spectrum of accelerated particles as a function of the velocity com- 
pression factor for different values of the shock speed. The scattering in the SPAS regime 
both upstream and downstream. 
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Fig. 5. — Angular distribution of accelerated particles at the shock front as measured in the 
downstream fluid frame, /i is the cosine of the angle between the particle direction and the 
shock normal. The different lines correspond to different values of the velocity compression 
factor, as indicated in the labels, while the shock velocity has a fixed value for each plot: 
^shPsh = 0.1 (top-left), 0.5 (top-right), 1.0 (bottom-left) and 5.0 (bottom-right). 



- 15 - 



the last paragraph of §3.21 A very similar behaviour was found also in (IKirk fc Duffy 19991 ) 



(compare their Fig 3 with left panel of our Fi g. El) where they use a different technique, 



the eigenvalue expansion firstly introduced by (IKirk fc Schneider 19871 ). and an equation of 



state for a gas consisting of both hydrogen and helium. 

Based on the discussion in §3.31 the effect of a turbulent magnetic field is expected to 
depend upon the shock speed. As illustrated in Fig. [7] (left and right panel as in Fig. [6]) 
the spectrum of accelerated particles is harder than in the absence of magnetic fields for 
newtonian and trans-relativistic shocks (the minimum slope is s = 3.75 in the equipartition 
regime, £ m ~ 1). In the ultra-relativistic regime all the curves approach the same value: the 
configuration of magnetic field adopted here does not produce any change in the particle 
spectra with respect to the case where no magnetic field is present. On the other hand this 
is easy to guess simply looking at Eq. 



The spectrum of the accelerated particles is however rather sensitive to the structure of 
the turbulent field in the downstream fluid. As we discussed in Sec. 13. 3[ the equation of state 
of the turbulent field depends on whether the field is generated downstream and is therefore 
isotropic in the local frame, or it is rather compressed in its perpendicular components. In 
this latter case the equation of state of the field is not the same as that of a relativistic fluid 
and this affects the compression factor at the shock. As an instance we consider here the case 
FshPsh = 5 and we introduce a magnetization parameter a = (5B 2 /8n)/ pc 2 in the upstream 
region, where 5B is the amplitude of the average magnetic turbulence, a = corresponds 
to the unmagnetized case, which leads to a compression factor 3.12 and a spectral slope 
s = 4.12, as already found earlier. On the other hand, for a = 10~ 2 (a = 3 x 10~ 2 ) the 
compression factor becomes 2.71 (2.16) and the spectral slope is s = 4.4 (s = 4.95). The 
corresponding value of the parameter £ m downstream, as resulting from the compression of 
the perpendicular components of the magnetic field is ^ m = 0.11 (£ m = 0.28). This softening 
of the spectrum may have very important phenomenological consequences for those classes 
of sources where particle acceleration occurs at ultra-relativistic shocks. 

Finally we consider the case in which both the effects of turbulent magnetic field down- 
stream and energy exchange between the thermal components of electrons and protons are 
taken into account. More specifically we concentrate on the so called equipartition case, in 
which £ e = £ TO = £ and we illustrate our results for different values of £. Here we restrict 
our attention to the case in which the turbulent field is generated downstream and does not 
result from the compression of an upstream field. As usual, the left panel in Fig. [8] refers 
to SPAS both upstream and downstream and the right panel to a large scale field upstream 
(and SPAS downstream). 

For £ = m e /m p the standard result is recovered in the SPAS case (left panel). The case 
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Fig. 6. — Slope of the spectrum of accelerated particles as a function of (3 s hT s h when the 
electrons have a temperature T e = £ e T p . The left panel refers to the case of SPAS in both 
the upstream and downstream plasmas. The right panel refers to the case in which a regular 
perpendicular field is present upstream. 




Fig. 7. — Slope of the spectrum of accelerated particles as a function of (3 s hT s h when a 
turbulent magnetic field is present with e m = £ m (e p — p p c 2 ). The left panel refers to the case 
of SPAS in both the upstream and downstream plasmas. The right panel refers to the case 
in which a regular perpendicular field is present upstream. 
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C, = when a large scale magnetic field is present upstream was discussed in (IMorlino et al. 20061 ) 
and the results are very close to the case £ = m e /m p used as the lower limit in Fig. [HJ 

In the other extreme case £ = 1 the spectrum of accelerated particles is harder than 4 in 
the non-relativistic limit and approaches the universal spectrum in the ultra-relativistic case. 
In both cases of SPAS and regular field upstream, a minimum in the slope of the spectrum 
is reached at f3 s hT s h ~ 0.4 corresponding to s ~ 3.65. This is a result of the interaction 
between electrons and protons in the trans-relativistic regime. 

In Figs. O [7] and [HJ the range of shock velocities adopted for the cases of upstream 
ordered magnetic field and of SPAS are quite different. This is not du e to any deep reason , 
but simply to the fact that the computation method introduced in (IBlasi fc Vietri 20051 ) 
becomes very challenging when the Lorentz factor becomes too large. This is related to the 
fact that in order to describe the regime of small pitch angle scattering we need to adopt a 
finite aperture a of the scattering function which howe yer must stay smalle r than 1/41^ in 
order to operate in the SPAS regime. As discussed in (IMorlino et al. 20061 ) this problem is 
absent in the case of ordered magnetic field. 



5. Conclusions 



We used the theoretical framework introduced by (IVietri 20031 ) and (IBlasi k, Vietri 20051 ) 
to determine the spectrum of particles accelerated at shocks with arbitrary speed and scat- 
tering properties, with different equations of state of the gas downstream. In particular we 
investigated two situations: 1) the downstream gas is made of thermal electrons and protons 
that may exchange energy with each other, thereby changing the equation of state; 2) the 
downstream gas includes a turbulent magnetic field. We also considered the case in which 
both effects are at work at the same time. 

In the downstream frame the scattering has always been assumed to be in the SPAS 
regime, while in the upstream fluid we considered two scenarios, namely SPAS and regular 
field with no turbulent scattering. We limited our attention to the case of strong shocks, 
namely the case in which the inflowing plasma has zero pressure. 

As a result of our calculations, we found several instances of violation of the so-called 
universality of the spectrum in both the newtonian and the trans-relativistic regimes. 

When the downstream plasma is made of electrons and protons and their temperatures 
are different, for newtonian and ultra-relativistic shocks the shape of the spectrum is not 
appreciably changed. However, for trans-relativistic shocks the compression factor increases 
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appreciably when the electron temperature is S> — T p , causing a flattening in the spectrum 
of the accelerated particles. In particular this is true for T e = T p , a situation which might be 
achieved due to some efficient collisionless plasma process, able to equilibrate the electron 
and proton components more efficiently than the simple isotropization of the velocity vectors. 

The effect of a turbulent magnetic field on the compression factor and on the spectral 
slope is more complex. We identified two situations of interest that arise when the magnetic 
field is introduced in the conservation equations at the shock surface: 1) the turbulent field 
is created downstream and is isotropic; 2) the magnetic field downstream is the result of 
the compression of the turbulent field upstream (only the perpendicular components are 
compressed). In the former case it can be demonstrated that the equation of state of the 
magnetic field is identical to that of an ideal relativistic gas, Eq. f[2"3l . irrespective of the 
shock speed. When the shock is non-relativistic the compression factor at the shock is 
increased and the spectra become harder than in the absence of field. For relativistic shocks 
the usual asymptotic spectrum is reached both in the case of upstream regular deflection or 
small pitch angle scattering. 

The second case is more interesting: when the perpendicular components of the tur- 
bulent upstream field are amplified by crossing the shock surface, the resulting magnetic 
field downstream is strongly anisotropic and the equation of state of the magnetic field is 
Pm = £m, not resembling that of a relativistic gas. In this case, for relativistic shocks the spec- 
tra of accelerated particles are softer than in the first scenario. For instance for T s h/3 s h = 5 
and a = 10~ 2 (er = 3 x 10~ 2 ) the compression factor becomes 2.71 (2.16) and the spectral 
slope is s = 4.4 (s = 4.95). The corresponding value of the parameter £ m downstream, 
as resulting from the compression of the perpendicular components of the magnetic field is 
£ m = 0-11 (Cm. = 0.28). Th e important effect consisting of a spectral steepening was found by 



(ILemoine fc Revenu 20061 ) and was attributed to the fact that the compression in the down- 
stream gas makes the magnetic field quasi-perpendicular, thereby reducing significantly the 
probability of return from the downstream frame. This effect is limited, in the analysis of 



(ILemoine &: Revenu 20061 ) to the particles with gyration radius smaller than the coherence 
scale of the turbulent field. The steepening of the spectrum as found in our calculations, is 
due to a change in the equation of state of the downstream plasma (electrons, protons and 
magnetic field) and concerns all of the spectrum of the accelerated particles. The relevance 
of this finding for the phenomenology of several astrophysical sources of accelerated particles, 
in particular those where a relativistic shock is expected or observed, is evident. 



The authors are grateful to an anonymous referee for the very interesting comments 
that helped us improving the manuscript considerably. This research was partially funded 
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Fig. 8. — Slope of the spectrum of accelerated particles as a function of (3 s hT sh with the full 
equation of state for the downstream plasma. The left panel refers to the case of SPAS in 
both the upstream and downstream plasmas. The right panel refers to the case in which a 
regular perpendicular field is present upstream. 



